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We report a numerical analysis of corrections to finite size scaling at the Anderson transition due 
to irrelevant scaling variables and non-linearities of the scaling variables. By taking proper account 
of these corrections, the universality of the critical exponent for the orthogonal universality class for 
three different distributions of the random potential is convincingly demonstrated. 
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The possibility of the Anderson localization of elec- 
tron states as a result of disorder was first suggested four 
decades ago H. Following the proposal of the scaling 
theory of localization M attention has focused on under- 
standing the critical properties of the Anderson transi- 
tion (AT), the quantum phase transition which occurs at 
a critical disorder separating a diffusive metallic phase 
from an insulating localized phase [gj. 

Our current understanding of the AT is based on the 
non-linear a model (NLcM) Q. This has been analyzed 
using an expansion in powers of e, where d = 2 + e is 
the dimension of the system. According to the NLcM it 
should be possible to classify the critical behavior using 
three universality classes: orthogonal, unitary and sym- 
plcctic depending on the symmetry of the Hamiltonian 
with respect to time reversal and spin rotation. Here 
we focus on the orthogonal universality class correspond- 
ing to systems with both time reversal and spin rotation 
symmetries. 

Beyond the suggestion of the appropriate universality 
classes, there has not been much success in making de- 
tailed predictions about the critical behavior with the 
NLerM. The problems are well illustrated by attempts to 
estimate the critical exponent v which describes the di- 
vergence of the correlation length £ at the AT. In early 
work it was found that v = 1/e |5j which gives v = 1 when 
extrapolated tod = e + 2 = 3. When combined with the 
Wegner scaling law s — ve J(| this leads to a conductiv- 
ity exponent s = 1. Measurements on some, but not all, 
materials do indeed yield s — 1 0. However, calcula- 
tions at higher orders in e produced strong corrections 
to the leading order when extrapolated to e = 1 H[§, 
showing that that this agreement is fortuitous. There is 
now no accepted estimate of the exponent based on the 
e expansion or any other analytic technique. 

While the above can be regarded as a rather unfortu- 
nate technical difficulty, fears have also been expressed 
that there maybe an infinite number of relevant opera- 
tors in the NLuM Jji| and that the theory may be un- 
sound. While it now seems unlikely in view of |11| that 
this is actually the case, in this context it is nevertheless 
important to have independent confirmation that our un- 
derstanding of the AT is correct. At present numerical 
simulations P,|l2]-|l4j offer the only viable alternative. 

In this paper we demonstrate an important basic prin- 
ciple underlying our understanding of the AT: the uni- 
versality of the critical properties of the AT. To do this 
has required us to address the principle uncertainty in 
previous numerical studies of the critical properties of 
the AT, the presence of systematic corrections to scaling 
in the numerical data due to the practical limitations on 
the sizes of system which can be studied. 

The computer time required in numerical studies of the 
AT increases very rapidly with increasing system size (as 
L 7 for the method used here.) This sets a severe limita- 
tion on the system sizes which can be simulated. How- 
ever, systematic corrections to scaling are expected in 
smaller systems and their neglect leaves important ques- 



tion marks over the validity of any conclusions drawn 
from the analysis of the numerical data. Here we con- 
sider two ways in which such corrections can arise: the 
presence of irrelevant scaling variables and non-linearity 
of the scaling variables [^5) . These effects lead to system- 
atic rather than random deviations from scaling and must 
be taken into account both when estimating the critical 
parameters and the likely accuracy of their estimation. 

Our work has also been inspired by the successful anal- 
yses of corrections to scaling in the Quantum Hall Effect 
(QHE) transition |IqjL7| ]. The present problem is, how- 
ever, more difficult since, unlike the QHE, the critical 
point is not known a priori on grounds of symmetry. 

The universality of the exponent for the box and Gaus- 
sian distributions of random potential was demonstrated 
to a limited extent in ]l2] | by taking account of corrections 
to scaling in an ad hoc manner. Here, taking account of 
corrections systematically, we confirm that result and ex- 
tend its validity to include the Lloyd model |l8j| . 

The Hamiltonian used in this study describes non- in- 
teracting electrons on a simple cubic lattice with nearest 
neighbor interactions only 
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Here i, y and z are the lattice basis vectors. The poten- 
tial V is independently and identically distributed with 
probability p(V)dV. We studied three models of the po- 
tential distribution: The box distribution 

p(V) = 1/W \V\ < W/2, 
= otherwise, 



the Gaussian distribution 
f 



p(V) = 



V2~7 



exp 



2cr 2 



with cr 2 = W 2 /12, and the Lloyd model in which V has 
a Lorentz distribution 
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ir{W 2 + V 2 )' 



For this distribution all moments higher than the mean 
are divergent and the parameter W is proportional to 
the full width at half maximum of the distribution. For 
these three models we analyzed the finite size scaling 
of the localization length A for electrons on a quasi-ld 
dimensional bar of cross section L x L. The length A 
was determined to within a specified accuracy using a 
standard transfer matrix technique [ ^9pc| ]. 

The starting point of our analysis is the renormaliza- 
tion group equation which expresses the dimensionlcss 
quantity A = X/L as function of the scaling variables 



A = / 
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In this equation b is the scale factor in the renormal- 
ization group, x the relevant scaling variable and ip the 
leading irrelevant scaling variable. We should find y < 
if ib is irrelevant. An appropriate choice of the factor b 
jL5| leads to 

K = F( X L 1 '\rbLy), (1) 

where F is a function related to /. 

For L finite there is no phase transition and F is a 
smooth function of its arguments. Assuming the irrele- 
vant scaling variable is not dangerous, we make a Taylor 
expansion up to order ni 

nj 

A = J2rL ny F n {xL^y (2) 

and obtain a series of functions F n . Each F n is then 
expanded as a Taylor series up to order ur 

F n { X L x ' u ) = J2 X m L m/ »F nm . (3) 
m=0 

To take account of non- linearities in the scaling variables 
we expand both in terms of the dimensionless disorder 
w = (W c — W)/W c where W c is the critical disorder 
separating the insulating (w < 0) and conducting phases 
(w > 0). 

XM = Er=l bnW n , = Enlo CnW". (4) 

The orders of the expansions are wir and to/ respec- 
tively. Notice that x(0) = 0. The absolute scales of the 
arguments in (Q) are undefined, we fix them by setting 
Fqx = F w = 1 in (^) . The total number of fitting param- 
eters is N p — (ni + l)(n,R + 1) + rriR + mj + 2 

The qualitative nature of the corrections can be under- 
stood by looking at some special cases. First let us sup- 
pose that non-linearities are absent (rriR — 1 and mj = 0) 
and truncate (^) at m = 1 

A = F (xL 1/v ) + $VFi (xL 1/v ) . 

From this equation we can infer that the estimate of the 
critical disorder, and possibly also the critical exponent, 
will appear to shift in a systematic way as the size of 
the system increases. To exhibit scaling it is necessary 
to subtract the corrections due to the irrelevant scaling 
variable. When nj = 1 we define 

Acorrcctcd = A - 1pL V F\ (x^) , (5) 

with the obvious generalization when nj > 1. We then 
have 

Acorrcctcd = -P± ( ~r ) ■ (6) 



The functions F± are defined by F±(x) = F (±(^±x) 1 ^). 
In this case the correlation length £ has a simple power 
law dependence on the dimensionless disorder £ — 
£± \ w \ V ■ The constants £± are not normally determined 
in finite size scaling studies. 

On the other hand, if we neglect the irrelevant variable 
and consider only non-linearity in the scaling variable 
we find A = F± (£/£) without the need to subtract any 
corrections. No systematic shift of the estimated criti- 
cal point should occur as the system size is increased. 
However, the correlation length £ no longer has a simple 
power law dependence on w but behaves as £ = £± \x\ 

The critical exponent v, the irrelevant exponent y and 
the functions F n are expected to be universal, while the 
coefficients {b n } and {c n } are not. Though we have ex- 
plicitly considered corrections due to the leading irrele- 
vant scaling variable only, the analysis can be easily ex- 
tended to several such variables. 

In the simulation A was evaluated as function of dis- 
order W for a range of system sizes L. The best fit was 
determined by minimizing the x 2 statistic pjj . This is 
justified if we suppose a uniform prior probability for all 
parameters, that the deviations between the model and 
the simulation data are purely random in origin and dis- 
tributed following a Gaussian distribution Q . This last 
assumption is also important in determining the likely 
accuracy to which the critical parameters have been es- 
timated. Since the inclusion of corrections to scaling al- 
lows for systematic rather than just random deviations 
from scaling in the numerical data, this assumption is 
more reasonable here than when corrections to scaling 
are neglected. Therefore we expect the estimates of the 
accuracy of the critical exponent etc. to be more reli- 
able. The model (|])-(|J) is nonlinear in some parameters 
so the goodness of fit Q has been checked using a Monte 
Carlo technique and the confidence intervals evaluated 
by re-sampling |22] , |2ljf| . 

The inclusion of the corrections in leads to a 

rapid increase in the number of fitting parameters and 
high quality data are essential if meaningful fits are to 
be obtained. All data used here have an accuracy of 
either 0.1% or 0.05%. To achieve this accuracy between 
10 6 and 10 7 iterations in the transfer matrix method were 
required. When deciding which correction terms to in- 
clude we attempted to maximize the goodness of fit Q 
while keeping the number of correction terms to a mini- 
mum. 

The details of the simulations and the types of fit used 
are listed in Table H. The estimated critical parameters 
and their confidence intervals are given in Table [n]. Some 
typical data are displayed in Fig. |l|. To exhibit scaling 
the data are re-plotted after subtraction of the appropri- 
ate corrections in Fig. ||. The corrected data now fall 
on a single curve clearly exhibiting scaling in agreement 
with (§). The magnitude of the corrections needed to 
obtain the scaling shown in Fig. || are of the order of 2% 
or so for the smallest system size decreasing to around 
0.3% for the largest system size. 
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The most important point to be drawn from Table |l| 
is that the estimates of the exponent v for the three dif- 
ferent disorder distributions are in almost perfect agree- 
ment. The same is true for the estimates of the critical 
parameter A c . This is strong evidence in favor of the 
universality of the critical exponent and other critical 
parameters. 

How do the results of the present analysis compare 
with those obtained when corrections to scaling are ne- 
glected? In Table [II we give the results obtained for the 
same potentials neglecting corrections. The first thing 
to notice is that range of system sizes (and in the box 
distribution, the range of W) for which an acceptable fit 
(Q > 0.1) can be achieved is very limited. After discard- 
ing data for the smaller system sizes reasonable agree- 
ment is obtained between Tables [n] and |lj|. However, 
given the more limited range of system sizes, the esti- 
mates of the accuracy to which the critical parameters 
have been determined when corrections are neglected are 
too optimistic. The problem is more evident when look- 
ing at less accurate e.g. 0.2% data 0| for the box distri- 
bution. Ignoring corrections to scaling, it was found that 
W c = 16.45 ±.01, A c = 0.586 ±.001 and v = 1.59 ±0.03. 
The estimates of W c and A c are not consistent with Ta- 
ble H The effect which gives rise to this inconsistency 
can also be seen in the data for Lloyd model displayed 
in Fig. [l]. A systematic shift of the apparent critical 
disorder to a lower value as the system size increases 
is evident. For the box and Gaussian distributions the 
shift was found to be in the opposite sense to higher dis- 
order. It seems likely that any analysis which assumes 
that deviations from scaling are purely random origin, 
rather than allowing for systematic corrections such as 
considered here, will lead to an over optimistic estimate 
of the accuracy to which the critical point has been de- 
termined and even to an incorrect determination of the 
critical point. In contrast the estimate of the critical 
exponent is quite consistent with that in Table [n]. Of 
course, the precise location of the critical point in any 
particular model is not in itself very important, but any 
inaccuracy in its estimate also affects the estimate of the 
critical conductance distribution p4| ]. 

We should also mention that surface effects and the 
influence of boundary conditions may also give rise to 
corrections to scaling behavior. We have used periodic 
boundary conditions to minimize surface effects. Even so 
there may remain some influence of the boundary condi- 
tions. To quantify this we have evaluated A for some 
representative values of the parameters for fixed (fbc), 
periodic (pbc) and anti- periodic boundary conditions 
(apbc). A large statistically significant shift in the local- 
ization length A was found between fbc and pbc. How- 
ever, even when calculating at a higher accuracy of 0.02% 
no such difference between pbc and apbc was found. We 
therefore think it reasonable to neglect such corrections 
here. 

We have presented a numerical study of the Anderson 
transition in three dimensions in which systematic correc- 



tions to scaling have been explicitly taken into account 
when estimating the critical disorder and other critical 
parameters. The universality of the critical exponent 
with respect to the choice of the distribution of disor- 
der has been accurately verified. 

While in this paper we concentrated on the scaling of 
the correlation length in the three dimensional Anderson 
model, corrections of a similar nature may also be im- 
portant in finite size scaling analyses of the conductance 
distribution and energy level statistics. The method we 
have described here is applicable in these cases and, in- 
deed, to any continuous quantum phase transition. 

Part of this work has been carried out on supercom- 
puter facilities at the Institute for Solid State Physics, 
University of Tokyo. 
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Disorder 


n R 


m 


m R 


mi 


w 


N d 


N p 


x 2 


Q 


Box 


3 


1 


2 





[15,18] 


224 


12 


214 


0.5 


Gaussian 


2 


1 


2 





[20.2,22.2] 


175 


10 


174 


0.3 


Lorentz 


2 


2 


1 





[4.1,4.5] 


224 


12 


203 


0.7 



TABLE I. The disorder distribution, the type of fit, the 
range of disorder W, the number of data Nd, the number of 
parameters N p , the value of x 2 for the best fit and goodness 
of fit Q. The system sizes used were L = 4,5, 6, 8, 10, 12, 14. 





Wc 


A c 


V 


V 


B 


16.54(53,56) 


0.576(74,78) 


1.57(55,59) 


-2.8(3.3,2.3) 


G 


21.29(28,31) 


0.576(74,77) 


1.58(55,61) 


-3.9(5.9,2.7) 


L 


4.265(52,72) 


0.579(76,88) 


1.58(47,65) 


-2.5(3.2,1.3) 



TABLE II. The best fit estimates of the critical disorder 
and the critical exponent and their 95% confidence intervals. 
The quantity A c = i*o(0) is expected to be universal. 





L 


W 


W c 


A c 


V 


B 
G 
L 


> 8 

> 8 

> 10 


[16,17] 

[20.2,22.2] 

[4.1,4.5] 


16.514(07, 22) 
21.28(26, 29) 
4.275(72, 78) 


0.579(78, 80) 
0.577(76, 78) 
0.574(73, 75) 


1.58(53,63) 
1.58(54, 62) 
1.58(53, 62) 



TABLE III. Best estimates of the critical parameters when 
corrections to scaling are neglected. 



FIG. 1. A as a function of disorder for the three dimen- 
sional Lloyd model. The solid lines are the fit (0)-([3). 



FIG. 2. The data in Fig. 1 after subtraction of corrections 
to scaling (see (JH|)) together with the the scaling functions 
(^). Here £ = £± |x| - "- The upper branch corresponds to the 
metallic phase and the lower branch to the insulating phase. 
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